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The deconfinement phase transition with external magnetic field is investigated in the Friedberg- 
Lee model. In the frame of functional renormalization group, we extend the often used potential 
expansion method for continuous phase transitions to the first-order phase transition in the model. 
By solving the flow equations we find that, the magnetic field displays a catalysis effect and it 
becomes more difficult to break through the confinement in hot and dense medium. 

PACS numbers: 24.85.+p, 25.75.Nq, 64.60.ae 


The most interesting quantum chromodynamics 
(QCD) phase transition at finite temperature and den¬ 
sity is the deconfinement from hadron gas to quark- 
gluon plasma (QGP). It is widely accepted that the 
quark matter may exist in compact stars with high 
baryon density and the initial stage of high energy nu¬ 
clear collisions with high temperature. Recently the 
study on the QCD phase transitions is extended to 
including magnetic field, since the strongest magnetic 
field in nature may be generated in relativistic heavy 
ion collisions. The initial magnitude of the field can 
reach eB ^ (1 — 100)m^ in collisions at the Relativistic 
Heavy Ion Collider and the Large Hadron Collider [ 11 - 
[^, where e is the electron charge and the pion 
mass in vacuum. When the magnetic field survives in 
the later formed hot medium, the interaction between 
the field and the medium may chaMe the fundamental 
QCD topological structure (see for reviews). 

Whether the magnetic field can induce deconfine¬ 
ment phase transition is still an open question. From 
the lattice simulations, the dependence of the Polyakov 
loop on the magnetic field supports the scenario of 
inverse magnetic catalysis, which shows a decreasing 
critical temperature as the field strength grows [M3- 
While in the MIT bag model and Polyakov quark meson 
model there is also the inverse magnetic catalysis m- 
[T3 - the magnetic field displays a catalysis effect in the 
Polyakov NJL model [l^] and turns to be an inverse 
catalysis only when a magnetic field dependent coupling 
constant is introduced [l3 • As for the phase transition 
at finite density, there is not yet precise lattice result 
due to the well-known sign problem. 

The Friedberg-Lee mode l Ha , also referred as a non- 
topological soliton model [^, is widely discussed in 
the study of confined and deconfined quarks [isl - l^ . 
In this model, the nonperturbative dynamics responsi¬ 
ble for confinement in QCD is simulated in terms of a 
nonlinear coupling of a scalar field cr. In vacuum state, 
the ensemble average of a is large and the quark mass is 
about 1 CeV, the heavy quarks are confined in hadron 
bags. With increasing temperature and/or density of 
the system, the average of a and in turn the effective 
quark mass drops down, the thermodynamic motion 
leads to the deconfinement of the effective light quarks. 

In this work we extend the Friedberg-Lee model to 
including an external magnetic field. We apply the 
functional renormalization group (FRC) method to the 
model to study deconfinement phase transition at fi¬ 
nite temperature, baryon density and magnetic field. 


As a nonperturbative method, the FRC [13113 is used 
to study phase transitions in various systems like cold 
atom gas [13113 1 nucleon gas [13 , hadron gas [13-113 > 
and quark matter [33 - 143 . By solving the flow equation 
which connects physics at different momentum scales, 
the FRC shows a great power to describe the phase 
transitions and the corresponding critical phenomena, 
which are normally difficult to be controlled in mean- 
field approximation because of the absence of quantum 
fluctuations. Instead of adding hot loops to the ther¬ 
modynamic potential in the usual ways of going beyond 
mean field, the FRC effective potential includes quan¬ 
tum fluctuations through mass, coupling constant and 
wave function renormalizations. 

The Friedberg-Lee model is defined as [13 

C = ip - ga)^!; + ^d'"adfj,a - U{a), 

U{cr) = ^ ^ ^ ^ ( 1 ) 


with ip, ijj, and a being quark, antiquark, and scalar 
fields, respectively, where we have neglected the cur¬ 
rent quark mass to simplify the calculations. There 
are four parameters in the Friedberg-Lee model, a with 
dimension L~^, b with dimension L~^, dimensionless 
c, and the Yukawa coupling constant g between quark 
and scalar fields. The bag constant B used to pro¬ 
vide a quark confinement is defined as the potential 
difference between perturbative vacuum and physical 
vacuum, B = C/((7per) — D((Tphy). The four parame¬ 
ters can be determined by solving the Euler-Lagrange 
equations and fitting the proton charge radius, proton 
magnetic moment and the ratio of axial-vector to vec¬ 
tor coupling [ilMlll. Different parameter sets are 
listed in Table U 221. 

Making a shift a = (cr) -I- 5a for the scalar field, 
where (cr) is the ensemble average, and substituting 
it into the Lagrangian density (Cl), the effective po¬ 
tential of the system at classical level is U{{a)), and 
the dynamical quark and sigma masses generated by 
the mean field (cr) can be extracted from the coef¬ 
ficients of the terms tpip and {5a)'^, nriq = g{a) and 
rn^ = a + h{a) +cl2{a)'^. The minimization of the clas¬ 
sical potential, dU{{a))/d{a) = 0, which is equivalent 
to the condition that the term linear in Sa disappears 
from the Lagrange density, is called the gap equation 
of the system and determines the value of (cr). From 
its solution, there are two vacuum states, the perturba¬ 
tive vacuum and physical vacuum, located respectively 
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TABLE I: The Friedberg-Lee model parameters [^ . 


set 

a(fm 


c 

9 

1 

69.73 

2112.6 

20000 

12.416 

2 

188.6 

7774.0 

100000 

15.333 

3 

11.60 

834.4 

10000 

10.957 

4 

17.70 

1457.4 

20000 

12.16 

5 

45.21 

5208.5 

100000 

16.379 

6 

6.734 

778.48 

10000 

10.963 

7 

10.25 

1358.4 

20000 

12.211 

8 

26.12 

4848.4 

100000 

16.537 


at (cr) Iper = 0 and (cr) |phy = ^ (l + - |p)- To 

simplify the notation, we take cr instead of (cr) in the 
following. 

We now apply the functional renormalization group 
method to the Friedberg-Lee model. The core quantity 
in the framework of FRG is the average action F k [cr, ip] 
The scale k dependence of the average action is charac¬ 
terized by the Wetterich flow equation [2^ [2^ in mo¬ 
mentum representation. 


With the chosen regulators Rl and i?^, after doing 
the three-momentum integration over the sigma and 
quark fields, the FRG flow equation ([2|) can be simpli¬ 
fied as 


dkUk 


tY— _ i _ 


-iN^NfTY 

n 


1 


( 5 ) 


where = 3 and Nf = 2 are color and flavor numbers 
of quarks, and Wm = 2m7rT and Wn = (2n -|- l)7rT with 
TO, n = 0,1, 2, • • ■ are the boson and fermion Matsubara 
frequencies in the imaginary time formalism of finite- 
temperature field theory. Doing the frequency summa¬ 
tion explicitly, the flow equation can be expressed in 
terms of the boson and fermion distribution functions 
nb{x) = 1/ — 1) and nf{x) = 1/ (e^ -|- 1), 

(1 - nf{Eqk + nf{Egk - /r)) , 

^ ^ak 


dk^k = :rTr 


-Tr 



-1 dRl 
dk 

cd 

+ 

[-H 

1_1 

dRl] 

dk 


( 2 ) 


where the trace Tr is defined in the inner and momen¬ 
tum space, and are the functional deriva¬ 
tives of the average action = S^Tk/Sa"^ and 

= STk!6{’ipip), and the infrared cutoff functions 
Rf. and Rl for scalar and quark fields, which are used 
to suppress quantum fluctuations at low momentum 
p < k, are chosen as the optimized regulators [4ll - l43| 


K = Y - P^) ® “ P^) 

^k = (7V+V7°) 


( 3 ) 


{po + ipY + k"^ 


(Po + ipf + 


- 10(fc2 -p2). 


Taking the local potential approximation, the aver¬ 
age action reads 


Ffe 


j d'^xi^ZliP - gko) ip + ^Z^d^ad^a 

+Uk{cj)), (4) 


where Z^, and Z^ are the wave functional renormal¬ 
ization constants, gk is the renormalized Yukawa cou¬ 
pling constant, and Uk{(j) corresponds to the potential 
in the Lagrangian density m but with scale dependent 
parameters Uk, bk, Ck and Bk- 

Assuming space-time independent ensemble aver¬ 
ages, and neglecting the Yukawa coupling constant 
renormalization and wave function renormalization 
{9k = g, Zl = Zl = 1), the average action to the low¬ 
est order is determined by the potential only [H, Issl - 
1 ^. IdSl - ld^ . Ffe = V 4 l 7 fc(cr), with the four dimensional 
space-time volume V 4 . 


where E^-k and Eqk are sigma and quark energies E^-k = 
yjk'^ + TO^ and Eqk = ^Jk!^ + with to^ = and 

TUq = ga. Given the initial potential Ga at the ultravi¬ 
olet momentum k = A, the potential at any k could be 
obtained by numerically solving the flow equation ([ 6 ]). 

In the Friedberg-Lee model, the quantum and ther¬ 
mal fluctuations come mainly from quarks which carry 
larger number of inner freedoms and lighter mass in 
comparison with the scalar meson. In the flow equa¬ 
tion ([ 6 ]) , the vacuum term and medium term of quarks 
have opposite signs, and the cancelation between them 
leads to the deconfinement phase transition. 

As a simple example, we consider the quantum fluc¬ 
tuations in vacuum where there is no thermal fluctu¬ 
ations. In the spirit of the Friedberg-Lee model, the 
sigma particle is heavy enough (rrin- > 2 GeV) to guar¬ 
antee the quark confinement [la, [i3]- In this case, its 
fluctuations in vacuum can be safely neglected and the 
flow equation becomes 


dkUk 


2 - — . 

TT^ Eqk 


( 7 ) 


Its analytic solution at fc = 0 can be obtained explicitly, 
I 


Uo{cr) = UA{cr) + 2AEqA (2A^ - 3(gcr)^) 


- 6 (gcr)'‘ln 


5CT 


A -t- EqA - 


( 8 ) 


with the initial condition [/a(o’). 

We now turn on the external magnetic field B in 
the Friedberg-Lee model and choose B = Be^ without 
loss of generality. Since cr is a neutral scalar field and 
there is no direct interaction with the magnetic field, 
we consider only the magnetic effect on the quark part 
of the flow equation. In the external magnetic field, 
the quark momentum integration is replaced by an in¬ 
tegration over the momentum along the z—axis plus a 
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summation over the Landau energy levels on the trans¬ 
verse plane 


/ 


d^p 





(9) 


with ai = 2 — 6io, electric changes Qu = 2e/3 and 
Qd = —e/3, and momentum p"^ = + 2\QfB\i. With 

this replacement, the flow equation dH) becomes 


dkUk = 


1 


1 


67r2 E^k V 2 

E EE 


27r^ 


s=±l f—u^d 2=0 


nb{E^k)j ( 10 ) 


E, 


qk 


0(fc^ -p\) (1 - nf{Eqk + /i) - nf{Eqk - p)) 


with quark transverse momentum = (2i -|- 1 — 
s)\QfB\ and the quark energy Eqk = ^Jk'^ + p\ Erriq. 

We solve the flow equation (nni) with the general¬ 
ized potential expansion method. The potential expan¬ 
sion method is successfully used in the study of chiral 
symmetry, C/jjT) symmetry and deconfinement at finite 
temperature |33l . l35l - l37l . l45l - l47l| . provided that the phase 
transition is continuous which guarantees only one min¬ 
imum of the potential. However, due to the character 
of the first-order phase transition in the Friedberg-Lee 
model, we should separately expand the effective poten¬ 
tial around the two local minima, and the comparison 
between them determines the real ground state of the 
system. This is a generalization of the potential expan¬ 
sion method from second- to first-order phase transi¬ 
tions. 

We expand the potential Ukicr) on the left-hand 
side of Eg. dTUl) around the two local minima U{ak) 
which satisfy the gap equation dUk((Jk)/dok = 0 and 
d^Uk{crk)/dal > 0. Shifting the field a ^ ak + Sak, 
dkUk could be expressed in terms of the powers of Sak, 


dkUkicr) 


1 dWk{<j) 

^ j\ da^ 
j=o 


{6akY ■ 

<7=ak 


( 11 ) 


By comparing the coefficients of {SakY on the left- 
and right-hand sides of the flow equation (uni, one ob¬ 
tains four coupled differential equations for the four 
parameters Ofc, bk, Ck and Bk- Note that in deriv¬ 
ing the four flow equations, we have used the relation 
dkSak = -dkCTk- 

Together with the gap equation, the k dependence of 
the four parameters and the classical field ak is fully 
determined. By solving the two sets of flow equations 
around the two local minima and comparing the ob¬ 
tained two effective potentials, the true ground state 
can be fixed. 

It should be mentioned that the flow equation for Bk 
should not be neglected when dealing with the first- 
order phase transition. According to the definition, the 
bag constant B is the potential difference between the 
two local minima. Since the two local minimal poten¬ 
tials are not fixed for a first-order phase transition, one 
should take into account its renormalization in the po¬ 
tential expansion. 


The initial condition for the four flow equa¬ 
tions and the corresponding gap equation at 
fixed temperature, chemical potential and mag¬ 
netic field is the four parameters and the clas¬ 
sical field at the ultraviolet momentum A, 
aA{T,p,B), &a(T,^, B), CA{T,fi,B), Ba{T,p,B) 
and <ta{T, p, B). Considering the fact that the system 
at high enough scale is dominated by the dynamics 
and not affected remarkably by the temperature, 
chemical potential and magnetic field, the tempera¬ 
ture, chemical potential and magnetic field dependence 
of the parameters at the ultraviolet momentum A 
can be safely neglected. Therefore, we take the 
initial values aA{T, p, B) = aA(0,0, 0), 6a(T,/ i, H) = 
6a(0, 0,0), ca{T,p,B) = ca(0,0, 0), Ba{T,p,B) = 
ilA(0,0, 0), and aA{T,p,B) = cta( 0,0,0). They are so 
chosen to reproduce the proton charge radius, proton 
magnetic moment and the ratio of axial-vector to 
vector coupling at fixed g in vacuum ai k = 0. The 
initial values corresponding to the different groups of 
physical parameters in Table [J are collected in Table [III 
at the ultraviolet scale A = 2000 MeV. For the results 
shown in the following, we take the first parameter set 
in Table H] with g = 12.416. For other parameter sets, 
similar conclusion can be obtained. 


TABLE II: The initial parameters in the FRG calculation. 


set 

aA(fm“^) 

1 

< 

1 

CA 

9 

1 

1749.51 

6400.35 

16549.2 

12.416 

2 

2733.81 

15263.9 

85174.6 

15.333 

3 

1307.08 

3457.19 

5148.48 

10.957 

4 

1610.3 

4950.51 

11692.3 

12.16 

5 

2928.46 

13477.2 

68802.3 

16.379 

6 

1302.34 

3369.6 

4805.52 

10.963 

7 

1614.94 

4854.62 

11117.1 

12.211 

8 

2963.68 

13287.9 

66526.8 

16.537 



FIG. 1: The effective potential at different momentum scale 
in vacuum. The dotted, dashed and solid lines correspond 
respectively to the scale k = 300, 268 and 0 MeV. 


How to choose the value of the ultraviolet scale A 
should be carefully discussed in effective models. In 
principle, the value of A should be large enough to guar¬ 
antee the saturation of the obtained physical result at 
k = 0. However, in models including hadrons as ele- 
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mentary constituents, the momentum scale can not go 
beyond the scale of the model itself where the hadrons 
are well defined. This means that the momentum scale 
should be restricted in a reasonable region. In the spirit 
of renormalization group, when more and more fluctu¬ 
ations are involved in the calculation through the mo¬ 
mentum scale k approaching from A to 0, the phase 
transition should happen at some critical scale kc- In 
the Friedberg-Lee model, this is shown in FiglT] At 
fc=300 MeV, the minimum of the potential is located 
at tr = a per = 0, the system is in the perturbative vac¬ 
uum state with deconfined quark mass mq = 0. This 
state maintains until k ^ kc = 268 MeV. The first- 
order phase transition happens at kc with two minima 
located at cr = aper = 0 and tr = aphy = 44 MeV. 
For k < kc, Uk<kr{^phy) becomes the true minimum 
and the system is in the physical vacuum state with 
confined quark mass niq = guphy = 547 MeV. The 
solid line is the final potential at fc = 0 with maximum 
fluctuations included in the calculation. From the dis¬ 
cussion on the phase transition in vacuum during the 
evolution of the flow equations, the ultraviolet momen¬ 
tum should be larger than the critical scale A > kc- 
Otherwise, there will be no phase transition at finite 
temperature, density and magnetic field. In the follow¬ 
ing calculation, we take A = 2000 MeV. 

FiglU shows the evolution of the dynamical quark 
mass mq as a function of scale k in vacuum. The 
quarks are massless before the critical scale k > kc, 
corresponding to the deconfined state. At the criti¬ 
cal scale kc, the mass jumps up suddenly, indicating a 
first-order phase transition between deconfinement and 
confinement. Note that, after the phase transition the 
quark mass is not a constant but changes gradually 
with the momentum scale k. This means that while 
the perturbative vacuum is always located at cr = 0, 
the location of the physical vacuum changes smoothly. 



FIG. 2: The quark mass rriq as a function of the momentum 
scale k in vacunm. 


Now we solve the flow equations at finite tempera¬ 
ture, density and magnetic field and determine the de¬ 
confinement phase transition point. Given fixed tem¬ 
perature, density and magnetic field, the dynamical 
quark mass is extracted from the evolution of the flow 
equations from fc = A to fc = 0. Figl3] shows the quark 
mass at fc = 0 as a function of the magnetic field at 
T = fj, = 0. With increasing magnetic field, the quark 


mass goes up continuously from 547 MeV in vacuum to 
606 MeV at eB = 16m^. The increasing quark mass 
indicates that, in an external magnetic field quarks are 
more tightly bound in hadrons and it becomes impossi¬ 
ble to break through the confinement by only magnetic 
field effect. 

Fig 0] displays the critical temperature Tc at vanish¬ 
ing baryon density and critical baryon chemical poten¬ 
tial He at vanishing temperature as functions of mag¬ 
netic field. They both increase monotonously with B, 
clearly indicating magnetic catalysis. 



eBIml 

FIG. 3: The quark mass niq as a function of magnetic field 
at zero temperature and baryon density T = fi ~ 0- 



FIG. 4: The critical temperatnre Tc at /r = 0 and critical 
chemical potential fic a-t T = 0 as fnnetions of magnetic 
field. 


In summary, the deconfinement phase transition 
under external magnetic field is investigated by apply¬ 
ing the functional renormalization group method to 
the Friedberg-Lee model. By expanding the effective 
potential around its two local minima and making 
comparison of them, which is a generalization of 
the usually used potential expansion method for 
continuous phase transitions, we determined the true 
minimum and the first-order phase transition point. 
Both the critical temperature and baryon density 
increase with the magnetic field strength, showing 
that the deconfinement phase transition in the model 
becomes more difficult in the external magnetic field. 
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